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A quantitative understanding of the dynamics of biological neural networks is fundamental to gaining 
insight into information processing in the brain. While techniques exist to measure spatial or temporal 
properties of these networks, it remains a significant challenge to resolve the neural dynamics with 
subcellular spatial resolution. In this work we consider a fundamentally new form of wide-field imaging for 
neuronal networks based on the nanoscale magnetic field sensing properties of optically active spins in a 
diamond substrate. We analyse the sensitivity of the system to the magnetic field generated by an axon 
transmembrane potential and confirm these predictions experimentally using electronically- generated 
neuron signals. By numerical simulation of the time dependent transmembrane potential of a 
morphologically reconstructed hippocampal CA1 pyramidal neuron, we show that the imaging system is 
capable of imaging planar neuron activity non-invasively at millisecond temporal resolution and micron 
spatial resolution over wide-fields. 



nformation processing in the brain is presumed to arise from interactions and correlations across several 
orders of magnitude of temporal and spatial scales and tens of thousands to billions of computational units. 
The smallest computational unit is the synapse, with sub -micron structures involved in chemical and electrical 
signalling. Changes in synaptic dynamics are the basis of learning and memory. While local conditions in the 
neuron determine signal flow into and out of synapses, complex signal integration and filtering occurs indepen- 
dently in different branches in the dendritic tree (the neuronal input structure). The decision to trigger an output 
is made in another structure - the axon. Neurons are organized in networks with complex, and largely unknown, 
connection rules. Networks may contain over a dozen different neuron types, each with their own dynamics and 
connection rules. Understanding the behaviour of these systems requires understanding the interactions between 
computational units at different scales in the system 1,2 . 

Current diagnostic techniques are limited in the number of computational units that can be recorded simul- 
taneously at sufficiently high spatial and temporal resolution. New techniques in optogenetics are enabling the 
manipulation of networks at some of these scales 3 however technology to read neuronal networks lags consid- 
erably. Voltage sensitive dyes allow readout of the neuron membrane potential, but have poor signal to noise 
properties and are toxic, making them unsuitable for long term recording 4 ' 5 . Voltage sensitive fluorescent proteins 
show promise, but also have poor signal to noise properties and limited temporal and spatial resolution 6 . Direct 
invasive techniques such as electrophysiological probes are fundamentally limited by the number of electrodes 
that can be placed in the tissue. Progress in silicon nanowire field-effect transistors have demonstrated sub 
millisecond temporal resolution, however spatial resolution may ultimately be limited by the distance between 
adjacent devices 7 " 10 . 

Our detection method is based on non-invasive magnetic field detection and is therefore fundamentally 
different to approaches based on detection of electric fields. While the techniques may ultimately complement 
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each other, our magnetic field based technique will not suffer from 
some of the drawbacks of electrically-based detection, such as probe 
positioning with respect to the Debye length and background electrical 
noise sources at the nanoscale. Furthermore, the detection is inher- 
ently non- invasive and issues such as the compatibility of quantum 
measurements of the NV defect with biological systems are now well 
established in terms of the low toxicity of diamond and the low photo 
and microwave powers involved. 

The detection set-up we consider consists of a commercial grade 
single crystal ultra-pure diamond membrane substrate containing a 
fabricated layer of negatively charged nitrogen -vacancy (NV) defect 
centres [Fig. 1]. The NV centre is a remarkable optical defect in 
diamond which allows discrimination of its magnetic sublevels 
through its fluorescence under illumination. Effectively, each NV is 
an atomic-sized magnetic field sensor that can be read-out either 
confocally 1112 or via a wide-field CCD 13 , and in a living cellular 
environment 14 . Neurons can be grown directly on the diamond sur- 
face 15 whose low toxicity is ideal for biological applications 16 . As we 
show, the ensemble of NV centres provides high sensitivity to the 
magnetic field fluctuations resulting directly from the transmem- 
brane potentials generated by the neural activity at sub -millisecond 
time-scales, and the spatial attenuation of the magnetic field at a 
100 nm standoff provides micron spatial resolution. This standoff 
is conservative, as implantation techniques permit the creation of NV 
centres within a few nm of the diamond surface. Employing a phys- 
ical model of the hippocampal CA1 pyramidal neuron, developed by 
Royeck et al 17 and modified by Wimmer et al 18 ' 19 , we show that the 
NV detection system is able to non-invasively capture the transmem- 
brane potential activity in a series of near real-time images, with 
spatial resolution at the level of the individual neural compartments. 
The data obtained will allow both the planar morphology and func- 
tion connectivity to be determined. The realisation of this detection 
system using available technology would represent a significant step 
forward in measuring and understanding the dynamics of whole- 
scale neuronal networks. 

In what follows, we first outline how transmembrane potentials 
generate magnetic fields, their typical strength and detection using 



the NV centre as a nanoscale magnetometer. We then analyse the 
detection sensitivity for a single axon case. We experimentally verify 
this sensitivity by propagating a simulated axon pulse along a micro - 
wire and detecting the resulting magnetic field signal with a prox- 
imate NV centre. Finally, we employ a model of a hippocampal CA1 
pyramidal neuron under a typical excitatory regime in which distal 
dendrites undergo current injections of roughly 2 nA. By direct 
simulation of the magnetic fields generated at the soma, axon and 
dendrites in response to this stimulation, we produce the corres- 
ponding image output (assuming current CMOS imaging techno- 
logy 20 ) and determine the effective spatial and temporal resolution of 
the system. 

Results 

Sensitivity Analysis. In order to establish our detection regime and 
required sensitivity we determine theoretically the magnetic fields 
generated by a transmembrane potential. We model an axon 
segment as a cylinder of radius a and length L aligned along the z 
axis [Fig. 2(a)]. These segments form the building blocks of more 
complicated neuron models to be considered later. Let €>(r, z, t) 
denote the radially- symmetric electric potential at radius r, 
longitudinal distance z and time t. Approximating the width of the 
cell membrane to be infmitesimally small, there is a step-change in 
the potential as it goes from just within the cell, 0(a_, z, t), to just 
outside the cell, €>(a+, z, t). The transmembrane potential, given by 
the difference V m (z, t) = <£(a_, z, t) — <£(a+, z, t), therefore 
represents the voltage drop across the cell membrane at 
longitudinal position z and time t. Using data for V m , one can 
reconstruct the 0(r, t) in the regions inside and outside a given 
component via the solution of Laplace's equation with boundary 
conditions set by V m (see Methods section). By solving for the 
potentials in both regions, we determine the electric field and 
hence current densities using Ohm's law, J(r,f) = <7(r,£)-V<D(r,£), 
which is integrated using the Biot-Savart law, 
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Figure 1 | Schematic of the nitrogen-vacancy (NV) in diamond neuron detection system showing a neural network component (in this case an axon) on 
a diamond substrate containing fluorescent NV centres, (a) The quantum state may be controlled via application of 2.88 GHz microwave radiation, and 
is dependent on the strength of the field produced by the axon, (b) Atomic lattice structure of the NV centre, (c). Upon optical excitation at 532 nm, the 
NV centre spin state may be measured (readout) by monitoring the intensity of the emitted red light using a CCD or CMOS camera, (d) Simulated 
dynamic output from a single CMOS pixel, (e) By monitoring many pixels, we may obtain a dynamic widefield image of the neural dynamics and network 
structure. 
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Figure 2 | (a)Schematic of the NV-diamond neuron detection system showing a neural network component (in this case an axon) on a diamond substrate 
containing fluorescent NV centres. For a FID based protocol, regions where the field is greatest produce the least fluorescence. For and ODMR based 
protocol, the converse is true, (b) Longitudinal variation of the transmembrane potential, V my of a crayfish lateral axon as taken from Ref. 21 . As this 
excitation propagates along the axon, it produces a time dependent magnetic field like that shown in (c) for a standoff of 100 nm. (d) Temporal 
arrangement of microwave and optical pulse sequences. For the FID based protocol, \ill microwave pulses are used to prepare the NV in a superposition 
of 10) and 1 1) states which accumulate a relative phase shift during the free evolution period due to the presence of a neural field. A second nil pulse 
transforms this phase into a population difference which is read out using a 532 nm laser. Application of this laser also acts to re-polarise the NV in the 1 0) 
state. The decreased intensity of the collected fluorescence permits the determination of the phase shift and hence the strength of the neural magnetic field. 
In the case of an ODMR based protocol, the constantly applied microwaves drive transitions between the paramagnetic sublevels of the ground states, 
whose relative energies are altered by the presence of a neural field, making the 10) — > 11) transition less likely, resulting in an increased fluorescence. 



to obtain the resulting magnetic fields, B(r, t). In Fig. 2 we show 
this for the case of a crayfish lateral axon using the measured 
transmembrane potential [Fig. 2(b)] taken from Ref. 21 . The 
resulting magnetic field signal calculated at the NV position is 
shown in Fig. 2(c) for a 100 nm standoff. Typically, the magnetic 
fields generated by the transmembrane potential are at the nT scale. 
We will show that the quantum detection system based on optically 
detected magnetic resonance (ODMR) of the NV centre in diamond 
has the combination of sensitivity and the appropriate spatial scale 
to enable detection of the neuronal magnetic fields at the high 
temporal-spatial resolution required. 

Theoretical proposals for the use of NV centres as sensitive, single 
spin magnetometers 22 ' 23 have been followed by proof- of-principle 
experiments 1112 ' 24 . Despite their promise, biological systems are gen- 
erally of a stochastic nature and are not always compatible with 
protocols designed to sense highly controllable DC and AC signals. 
It was recently shown that decoherence microscopy methods 25 permit 
the characterisation of randomly fluctuating magnetic fields 26 . 
Methods to characterise random AC fields of random frequency 
and phase have also been proposed 27 . In addition to their use as 
fluorescent biomarkers 16 ' 28 " 31 , NV centres have recently received 
much attention as sensors of magnetic biological processes, in par- 
ticular as sensors of ionic flux through single ion channels in lipid 
membranes 32 , and intra- cellular micro-fluid dynamics 14 . The advent 
of widefield NV magnetometry based on CCD detection of a high NV 
density diamond substrate has yielded the ability to spatially recon- 
struct magnetic fields at the fim scale together with improved sens- 
itivity over that obtainable with a single NV centre whilst retaining 
the necessary quantum coherence properties 33 . Further improve- 
ments are expected in the near future, with nitrogen to NV~ conver- 
sion efficiencies in excess of 50% having been recently reported 34 ' 35 . 

We envisage two detection protocols based on either directly mon- 
itoring the coherent phase difference between the sublevels of the NV 
centre using free induction decay (FID) experiment; or by measuring 
the location of the resonance peak in a continuous wave ODMR 



experiment [Fig. 2(d)]. The sensitivities of the two methods are 
essentially equivalent 24 , however for defmiteness (and brevity) we 
describe in more detail the former. At fiducial t = 0, following optical 
polarisation into the m s = 0 state, the application of a nil microwave 
pulse places each NV centre into an even superposition of m s = 0 and 
m s = 1 Zeeman levels. After a free evolution time of t, another nil 
pulse is applied and the system's state is read out optically. In the 
absence of any magnetic disturbances the system remains coherent, 
and the probability of finding the NV system in the m s = 0 magnetic 
level after the second nil pulse, measured directly from its fluor- 
escence, will oscillate between 0 and 1 at a frequency given by the 
energy difference of the two levels. However, in reality such coherent 
superpositions of quantum states are highly sensitive to magnetic 
field fluctuations in the immediate crystalline environment due to 
mutual spin flipping of electronic nitrogen defects and/or 13 C nuclear 
spins. Thus, depending on the material composition, one typically 
measures a decay of these coherent oscillations between the magnetic 
levels as a function of the evolution time t. For isotopically pure 
diamond the timescale of this decoherence, referred to as the FID 
time TJ, is typically of order 1 to 10 fis and ultimately sets the 
sensitivity limits for an NV based magnetometer. By employing a 
spin -echo pulse sequence in which an additional n -pulse is applied at 
time t/2, coherence times may be extended by more than 2 orders of 
magnitude, and further improvements may be realised by employing 
higher order decoupling schemes such as CPMG 23 ' 36 or Unrig 37,38 
pulse sequences. Such decoupling schemes achieve improvements 
to the sensitivity by suppressing the effect of low frequency noise 
on the NV centre, and are hence not suited to measuring neural 
magnetic fields which fluctuate on T = 1 to 10 ms timescales. To 
ensure maximal sensitivity to low frequency noise we must employ 
some N t repetitions of the FID or ODMR protocols in the time 
interval [£, t + St]. As such, St serves to define the temporal resolu- 
tion of the system, and must satisfy St < T. Furthermore, each 
protocal cycle must be sufficiently sensitive that N t < T /T^ mea- 
surements may be taken within the timescales associated with the 
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neuronal dynamics, whilst maintaining an acceptable signal to noise 
ratio (SNR). 

During its evolution under the neural magnetic field B(r, t), the 
coherent phase accumulation is of an NV centre at position r, during 
time interval t, is given by A0(r,-) = y J Q T Bi(r i9 t)dt, where y is the NV 
gyromagnetic ratio and Bj is the component of B parallel to the i th N- 
V axis. The signal change we are interested in measuring is the 
resulting change in fluorescence due to this phase accumulation, 
which is proportional to sin (A0). Since the neural magnetic fields 
considered here are typically less than 10 nT, and are essentially 
constant over the timescale of a single measurement protocol, we 
have that A0<O, thus the singal change at a single defect site is given 
by Si~yTBi(Yi)V(Yi,T), where V(r^%) is the decoherence envelope of 
the NV defect at position r f as a function of the interrogation time, t. 
This envelope represents the decay of quantum coherence of defect i 
due to interactions with the surrounding magnetic environment, 
including 13 C nuclei, nitrogen electron spins, and other proximal 
NV centres. As we are considering an FID sequence, the shape of 
each envelope is assumed to be Gaussian, and is given by 



V{t) = exp 



(rV) 2 -(r>) : 



(2) 



where P = T (r f ) is the intrinsic decoherence rate of NV centre i, and 
r Cl3 ~ 100 kHz. The decoherence rate due to nitrogen impurities, T n 
will explicitly depend on the nitrogen density and must be optimised, 
as shown below. The total signal from a single FID protocol cycle, 
S= Si, is the sum of the signals from all NV centres, and in the 
high NV density limit, is given by 



S(t) = dyxn Bi(r,T)V(r,T)d 3 r, 



(3) 



where n is the total electron spin density due to both nitrogen and NV 
defect centres (assumed uniform for chemical vapor deposition 
grown samples, and an axially- symmetric Gaussian for implanted 
samples), a is the nitrogen-NV conversion efficiency, and V p is the 
volume of a given pixel. Since the measurement of each NV centres is 
projective, the total signal distribution follows the central limit the- 
orem and the uncertainty in the signal is given by 



3S-- 



(4) 



1 

C\ N t 9 

where C is an experimental parameter accounting for imperfect 
photon collection (including non-unity quantum efficiency of the 
detector) and signal contrast, and N p is the effective number of NV 
defects (probes) in V p , given byN p ~ anV p . To measure the magnetic 
field of an axon, ultimately we require S > 3S. 

We now consider a simple, analytically solvable case where the 
axon dimensions are large (~ 60 ^m, as in the case of a crayfish 
lateral axon 21 ) compared with a sensing volume of V p ~ (1 fi m) 3 , 
ensuring there is little variation in the magnetic field strength over 
the sensing volume (i.e. B { = BVi). The signal is then S ~ OLnV p yBx. 
The minimum detectable magnetic field is then 



yC^LNtV v 



ty/n 



which we wish to optimise for n and t, giving firstly 

/ ?\ -1 / 2 
?opt = ( + 2T^ J . The dephasing due to nitrogen is given 



by r n 



Kn 38 , where k - 



~hf. Setting £SB = 0 gives 



r ci3 =r n , hence n = T Cl3 /K = 23 x 10 23 m -3 = 1.3 ppm, an order 
of magnitude below the maximum reported density achieved in 
practice to date of 2.8 X10 24 m" 3 39 . The minimum detectable field 
by a pixel of volume V p is then 



3B = 



3.3 nT ^m 3 / 2 



For example, an implementation with C = 0.3 and an N-NV con- 
version ratio of a = 0.1 would be capable of resolving a magnetic field 
to an accuracy of 3.5 nT in an integration time of 1 ms (N t = 100). 

In the large axon case considered here, increasing V p will result in 
an improved sensitivity, however this will not be true in general. As 
V p becomes comparable to the axon dimensions, the field felt by 
distant NV centres will be significantly less than that felt by those 
proximate to the axon, and the integral in Eq. 3 will no longer scale 
linearly with V p . The noise amplitude however [Eq.4] will grow with 
the square root of the sensing volume, regardless of the axon field 
characteristics. If the sensing volume is sufficiently large that the 
signal exhibits sub -square root scaling with the sensing volume, there 
will be no advantage in having a larger V p . As such, better results for 
smaller neural components will be achieved by optimising the V p for 
the task at hand. 

The maximum frequency with which measurements maybe taken 
is fm = (T opt + 0 _1 > where T opt = (2r c ) _1 from above, and T m ~ 
350 ns is the time required for photon collection and subsequent re- 
polarisation of the NV spin state. To gain further improvements to 
the SNR, we envisage taking numerous measurements with fre- 
quency / m , and applying a low-pass filter to the measurement record. 
Incorporating a greater number of timepoints in this process will 
improve the sensitivity, however this comes at a cost of decreased 
temporal resolution, 3t = N t /f m . A faithful reconstruction of the field 
dynamics with time-lag 3t will be possible provided 3t is less than the 
characteristic timescales of the neural dynamics. 

2 



dt = 



3.3 nTV 



L opt 



KVp/mi 3 / 2 



(6) 



The effect of averaging over different ranges of measurements and 
the trade-off that exists between a high SNR and high temporal and 
spatial resolution is considered later for the case of the Hippocampal 
CA1 pyramidal neuron. 

Comparison of magnetic and electric field detection. This method 
of monitoring neuron function through the magnetic fields 
generated offers a distinct advantage over consideration of the 
resulting electric fields. The electric field is a consequence of the 
local gradient of the electric potential, whereas magnetic field 
detection is sensitive to non-local field sources. That is, the 
external magnetic field is a consequence of both the internal and 
external electric fields, the strength of the former being some 3 
orders of magnitude larger than the latter. 

Electric field detection protocols using single electron transistors 
(SETs) have sensitivities of 2 Vein" 1 Hz" 1/2 at standoffs of 100 nm 41 , 
however there are significant practical compatibility issues associated 
with real neural samples and the need for cryogenic sensing opera- 
tion. Electric field sensing using NV centres does allow for room 
temperature operation, however the associated DC field sensitivity 
was recently demonstrated at 613 V cm" 1 Hz" 1/2 42 , whereas single 
NV sensitivities to DC magnetic fields have been demonstrated at 
43 nT Hz" 172 40 . The peak magnitudes of magnetic and electric fields 
around a typical axon are shown in Fig. 3(a), yielding 3 nT and 
40 X 10" 3 V cm" 1 , respectively, at a standoff of 100 nm. For the sake 
of comparison, Fig. 3(b) shows the time required for a single NV 
centre to resolve the magnetic and electric field strengths in the 
region surrounding a typical axon. Clearly local magnetic detection 
offers an advantage in terms of both spatial and temporal resolution. 

Experimental verification of detection protocol. To experimentally 
demonstrate the effectiveness of the proposed detection protocol for 
neuron specific signals, and quantitatively verify the theoretical 
analysis, we replicated the magnetic field produced by a single 



SCIENTIFIC REPORTS | 2 : 401 | DOI: 1 0.1 038/srep00401 



4 




Radial distance above axon surface [m] 



Figure 3 | (a) Dependence of magnetic and electric field strength on radial distance above the axon surface. The more immediate attenuation of the 
magnetic field permits roughly an order of magnitude better spatial imaging resolution as compared with the electric field. Note the different axes used for 
electric and magnetic fields, (b) Comparison of the temporal resolution corresponding to sensing the magnetic field using diamond crystal containing a 
high NV centre density (NV, blue), and sensing the electric field using the same NV centre density (NVE, green) and a single electron transistor (SET, red). 



axon using a microwire on the surface of a diamond substrate 
[Fig. 4(a)]. The pulses were constructed such that the resulting 
Biot-Savart field emulates the temporal dynamics of the axon 
considered in the preceding section [Fig. 4(b)]. A single NV centre, 
located at a lateral distance of approximately 10 /urn from the wire, 
was used as the sensor (see Methods section). 

We simulated the detection gain of N p centres by repeating the 
measurement process N s times. The binning widths for the photon 
integration times were t p = 24.2 fis, and the total time taken to 
send two successive pulses down the wire was 20 ms [Fig 4(d)]. 
The dynamics of this signal place a limit on the permis sable 
temporal resolution, St ~ 1 - 2 ms, hence the maximum number 
of data points employed in the low-pass filtering process is N t = 
dtltp ~ 50. Assuming a constant sensitivity, the effective number 
of centres involved in the detection may then be determined from 



N p = — ( — ] , where B w is the peak magnetic field strength due 

N t \B a J 

to the pulse in the wire, and B a is the peak magnetic field strength 
due to the axon. This assumption was experimentally verified by 
measuring the minimum detectable magnetic field for a range of 
photon integration times [Fig 4(c)], from which a corresponding 
sensitivity of 10 f/T Hz" 172 was determined [Fig 4(d)], for this 
(non-optimised) system. 

The figures of merit are the effective spatial and temporal resolu- 
tions, given by Sx = (N p /n) l/3 and St = N t t p respectively. The mea- 
surement record obtained by monitoring a pulse with B w = 24 f/I, 
and N s = 5, 000 cycles is given by the green trace in Fig 4(b). From the 
above scaling, this is equivalent to a pixel volume of V p = (7.9 /mi) 3 at 
a temporal resolution of St = 1.3 ms. Sacrifices in the temporal 
resolution will allow for decreases in the required pixel volume, for 
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Figure 4 | Detection of experimentally simulated neuron signals using the ODMR detection protocol on a single NV centre, (a) Confocal fluorescence 
image showing the spatial arrangement of the NV with respect to the microwave and artificial neuron signal micro-wires, (b) Experimental measurement 
of the electronically generated neuron pulse together with the extrapolated effective spatial and temporal resolution corresponding to a high NV density 
widefield detection system. Optimisation of this protocol to achieve higher resolutions is discussed in the main text, (c) Experimental minimum 
detectable magnetic field, SB vs total photon integration time, T, showing a SB oc T~ 1/2 dependence. This dependence is verified by plotting the 
corresponding sensitivity, rj = SBy/T, in (d), showing the sensitivity is effectively constant. 
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example, a pixel volume of V p = (1 /im) 3 has a corresponding tem- 
poral resolution of St = 2 ms. These figures clearly demonstrate, 
experimentally, that NV are centres are capable of simultaneously 
resolving both the spatial structure and temporal dynamics of neur- 
onal magnetic fields, even in their current, unoptomised imple- 
mentation. Further improvements may be achieved with the use of 
higher grades of isotopocially pure diamond crystal, as used in 40 , 
where single-spin DC magnetic field sensitivities of 43 nT Hz" 172 
were demonstrated. Extrapolating this to the present context, such 
sensitivities would permit a pixel volume of V p = (0.2 fim) 3 at a 
temporal resolution of 1.3 ms. How these capabilities relate to the 
spatial and temporal characteristics of a biological neuronal network 
are discussed in the following section. 

Imaging simulation: the hippocampal CA1 pyramidal neuron. In 

order to quantitatively describe how the device would sense and 
image neural activity we simulated the magnetic fields that would 
be produced by a neuron while it receives synaptic input and 
generates and fires an action potential output. The neuron model 
we used was of a morphologically reconstructed hippocampal CA1 
pyramidal neuron. Hippocampal CA1 pyramidal neurons have an 
extended morphology, show a rich repertoire of dynamics and are 
involved in networks that underlie important behaviour such as 
learning and memory. Study of these neurons and the networks 
they reside in will be a major target of the techniques described here. 



Typically such an experiment would be performed in vitro using 
a buffer solution of phosphate buffered saline (PBS), or equivalent. 
Such solutions contain dissolved sodium and potassium salts at 
concentrations of roughly 150 mmol L" 1 , giving rise to nuclear 
spin concentrations of roughly 10 25 m" 3 , some 5 orders of mag- 
nitude less than the hydrogen nuclear spin concentration due to 
water molecules. The hydrogen spins themselves produce an 
effective stochastic magnetic field with fluctuations caused by 
molecular self-diffusion. This would cause the decoherence rates 
of NV centres at even a few nanometres from the diamond surface 
to increase by approximately 100 Hz-1 kHz 32 , some 2-3 orders of 
magnitude less than that due to intrinsic decoherence sources, and 
can therefore be ignored. 

The model consists of 265 anatomical sections [Fig. 5(a)] and 15 
voltage and calcium activated conductances with a non-uniform 
distribution across the morphology. We stimulated the model neu- 
ron with 2 nA current injections into 21 sites on distal dendrites. The 
associated transmembrane potentials at specific locations are shown 
in Fig. 5(b). As the distal dendrites are stimulated current flows along 
these processes towards the soma. The associated magnetic field 
strength changes are readily seen in the sequence of time snapshots, 
Fig. 5(c-e). As the membrane potential increases the soma current 
flow along the dendrites decreases, reducing the magnetic field 
strength [Fig. 5(c)]. As active conductances in the soma and axon 
are recruited, current flow increases in these compartments and the 
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Figure 5 | (a) Morphology of the Royeck-Wimmer model and input excitation sites, (b) Transmembrane potentials generated at the soma, axon and 
dendrites shown in (a), (c-e) Zoomed plots of magnetic field strength at 100 nm standoff showing the integrate and fire effect of the central soma and the 
reactionary dynamics in the dendritic region below. Simulated measurements taken at (f) 0 /urn and (g) 45 /urn along the apical dendrite directly above the 
soma, as shown in (a) for a detection volume of V p = (2 /mi) 3 and a range of integration times, St. Note the change of magnetic field strength scale in (f) 
and (g), resulting in the latter requiring longer integration times. Assumed parameter values are C = 0.05, and a = 0.1. 
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field strength increases. Eventually an action potential is triggered in 
the soma and distal axon [Fig. 5(d)] . The initiation point in the distal 
axon is not as bright as the soma because the technique is sensitive to 
longitudinal currents, not voltage. A back propagating action poten- 
tial is also triggered in the basal dendritic tree [Fig. 5(e)]. 

In Fig. 5(f)&(g) we show the single (2 fim) pixel detection trace of 
the neuronal magnetic field at two locations below the neuron 
structure in response to the applied stimuli, assuming experimentally 
realised values of C = 0.05 43 , and a = 0.1. The magnetic field signal is 
plotted together with the simulated measurement output for a range 
of integration times, St. As the integration time is increased the SNR 
improves, at the cost of temporal resolution in the neuronal signal 
itself. However, it is evident that a minimal time-lag of 1 ms gives an 
acceptable account of the neuron signal in all three cases. 

Finally, we determine the overall spatial-temporal resolution of 
the imaging system by explicitly considering the trade-off 
between integration time and detection volume. In Fig. 6 we show 
the simulated CMOS/CCD output (assuming a readout rate of 
500 fps, as available with current technology 20 ) for a snap-shot of 
the neuron activity at t = 157 ms for a range of integration times 
St and spatial detection volumes Sx. It is clear that the neuronal 
magnetic field structure is apparent at [Sx, St] = [2 /urn, 1 ms], 
showing the imaging system has the temporal resolution to fully 
map the magnetic field dynamics of a neuronal network, whilst 



simultaneously reproducing the structural morphology at the 
sub -cellular level. 

Discussion 

We have investigated the use of the magnetic field sensing properties 
of NV-ensembles in diamond for imaging neural activity. Using 
published crayfish lateral axon biophysical data we determined the 
magnetic field signal at 100 nm from the axon surface to 
demonstrate that it lies within the sensitivity range of our detection 
system. The sensitivity regime for detection of the magnetic field 
dynamics generated in the axon structure was determined based 
on both pulsed (FID) and continuous wave (ODMR) sequences. 
Direct measurement of axon -scale magnetic fields, produced by pas- 
sing current in a proximate microwire, verified that the sensitivity 
limits fell within the range needed to detect neuronal signals. To 
explore the potential utility of NV arrays as wide field detectors of 
neuronal network activity we simulated the three dimensional mag- 
netic fields associated with action potential propagation in a mor- 
phologically realistic hippocampal CA1 pyramidal neuron placed 
100 nm from the NV detection layer. The simulated photon emis- 
sion of our model neuron/NV-layer combination was projected onto 
a virtual CMOS array to investigate the performance of the sensor 
in detecting local and wide-field neuronal structure and dynamics. 
By exploring different combinations of integration times, St, and 



St = 1ms St = 4ms St = 10ms 




Figure 6 | Simulated snapshot of the CCD output at t = 157 ms and for a range of detection volumes, V p = (<£x) 3 and integration times St. Assumed 
parameter values are C = 0.05, and a = 0.1. 
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detection regions, Sx, we found that the performance of the sensor 
enabled simultaneously high temporal and spatial resolution of 
extracellular field potentials in regimes beyond those obtained by 
current methodologies. In summary, our experimental results and 
theoretical work establish the significant potential of this quantum 
based technique to visualise the key components of neuronal net- 
work activity, subthreshold signalling, action potential initiation and 
propagation in axons, soma and dendritic compartments, at relevant 
scales to provide new views into network function. To bring this 
imaging concept to reality one must assemble all the individual 
and non-trivial components, such as sufficiently dense near-surface 
NV layer structures in high grade diamond material, neuron growth 
on these diamond surfaces, and microwave quantum control and 
optical wide-field readout commensurate with long term integrity 
and function of the biological structure. Precursor experiments based 
on well characterized biological systems will extend the experimental 
work carried out here on model neuronal signals. 

Methods 

Optically detected magnetic resonance (ODMR). ODMR detection is based on 
monitoring the NV fluorescence spectrum in response to an applied oscillating 
microwave field. The constantly applied microwaves drive transitions between the 
paramagnetic sublevels of the ground states, whose relative energies are altered by the 
presence of a background field. The effect of the field, B, will be to shift the frequency 
spectrum, S(co), by an amount yB to give S(co) =S(co — yB). At some chosen reference 
point in the spectrum there will be a corresponding change in the number of photon 

counts proportional to AN (go) oc AS(eo) = S(co) — S(ooJ. For small fields, this change will 

depend on the local gradient of the spectrum, AS— — S(co) Aco = yB — S(co) — KyB. 

aco aco 

The spectrum is given by a Lorentzian distribution of Full Width Half Maximum T, 

centred about yB 0 , 



S(co)=N 0 [ 1 



r 2 +y 2 (B-B 0 y 



(7) 



where N 0 is the unperturbed photon count in regions of the spectrum far detuned 
from resonance and C, is the difference in photon counts between the resonant and 
unperturbed regions. The maximum gradient occurs at go = +yJBo ± T/ V?>, giving 



= max< —S(co) 
I aco 



= CN 



3^3 
~8T~' 



(8) 



To ensure a signal to noise ratio of greater that 1, we require AS > \j \fN =\j \/Rt, 
where N is the total number of photons collected in time t, and R is the photon count 
rate. Hence the minimum detectable field strength and corresponding sensitivity are 

" 1 8r m 



CyVRi3V3' 

1 sr 

CyVR3V3' 



(10) 



The countrate, R, is a monotonically increasing function of the laser power that satu- 
rates at approximately 140 kHz, and both C and T are functions of the microwave and 
laser power. Maximal sensitivity will be achieved by optimising these parameters to 
reach the operating regime of the ODMR protocol that yields the maximum value of k. 

Calculation of the magnetic field from a transmembrane potential. Numerical 
simulations of the transmembrane potentials of the Royeck model were performed 
using the NEURON software package. Given an arbitrary form for the 
transmembrane potential along a given component of the neural network, we may use 
the corresponding Fourier transform, J 7 [V m (z,t)], with respect to the longitudinal 
coordinate z and conjugate variable k to compute the electric potentials inside and 
outside of the component. 

- r[Vm M\kwm«)~\ 



<S> t (r,z,f)- 



2n 



G(\k\a) 
M\k\r)h(\k\a) 



G(\k\a) 



where o { and o Q are the conductivities of the interior and exterior regions respectively, 
and I n and K n are the n th order modified Bessel functions of the first and second kind 
respectively, T~ x denotes the inverse Fourier transform, and G = GfK\h + ^i^oh- 
Current densities may be found from Ohm's law, f ° = o it 0 VO' °. For the cases 
considered here, the longitudinal current density is approximately 2 orders of 
magnitude larger than that in the radial direction, which we consequently ignore. The 
longitudinal current densities are given by 



kT\V m _ 



G(\k\a) . 
,^o(Wr)/i(Wa 



G(\k\a) 



The magnetic fields due to the current densities may be evaluated using the Biot- 
Savart law, 

Ho fj(r')x(r-r') , 3 j 



Kr)-^ f J(r/ ) x(r ~ 
U 4ti J | r -rf 

x {\- X 2 )K(x 2 )-2E(f) ^ 



where B^ is the azimuthal component of the field, and K (% 2 ) and E (% 2 ) are the 
elliptical integrals of the first and second kind respectively given by 



K(x 2 )= f"[l-Z 2 sin 2 (a)] 
Jo 

E(X 2 )= f"[l-f sin 2 («)] 1/2 do< 
Jo 

2 4ar 

and 7 = ^ 7 , 

{a + rf + {z-z'f 



Experimental verification of ODMR based measurement protocol. Following 21 , we 
fit the transmembrane potential using a sum of three Gaussian curves, 

V m = V 0 ^2 ^-exp [ - B] (f - df] , 

hence the external magnetic field near the membrane will be proportional to 

^ = -21/ 0 ^B 2 ( t -C 1 )A,exp[-B 2 ( t -C,) 2 ]. 

The experiment was performed by using an arbitrary waveform generator to 
propagate two sequential pulses of this shape along a 10 /mi wire. The resulting 
magnetic field was measured by the change in fluorescence of an NV centre separated 
from the wire by a lateral distance of approximately 10 /mi. The dependence of the 
magnetic field strength on the signal voltage was determined by driving DC signals in 
the wire and monitoring the peak splitting of the NV ODMR spectrum. From this, we 
determined the NV centre experiences a magnetic field of 120 nT for every mV 
applied to the wire. This corresponds to 360 nT for every mA of current in the wire 
due to an electrical resistance of 3 Q. 

Since the magnetic field is proportional to the applied voltage, we use the following 
fit for the NV fluorescence 



S = S 0 - 



;=1 i=l 



x exp 



(ii) 



where S ; , tj and t ; - are the relative amplitude, offset and dilation of the/ h pulse, and A t 
Bj and Q are as above. An example of the fit to a pulse with peak voltage of 200 mV is 
shown in Fig. 4(b). Each timepoint represents a photon integration time of 24.2 /is. 
The effect of increasing the number of NV centres involved in the detection is 
simulated by varying the number of cycles per time point from 1,000 to 50,000. 

From this fit, we may determine the mean square error, which gives the uncertainty 
in the photon counts, SN. The maximum signal contrast, AN is given by the difference 
between the maximum value of the fitted function, max(S) and the value that 
represents 0 magnetic field, S 0 . That is, AN = max(S) — S 0 . The uncertainty in the 
magnetic field strength is defined such that the magnetic field measurement has the 
same signal to noise ratio as that of the photon counts, 

AB AN 

Since we know the peak magnetic field is 120 nT per mV of signal in the micrwire, we 
may determine the corresponding uncertainty in the magnetic field, which in turn 
defines the minimum detectable magnetic field, 



5B = AB 



3N 
'AN 



--AB- 



SN 



max(S) — S 0 
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Summary of physical and detection timescales. 



Typical neronal pulse duration 
NV spin dephasing time 
Peak neuronal magnetic field value at 
100 nm from axon surface 
Single NV centre sensitivity 

High density NV centre sensitivity 
Maximum CMOS/CCD framerate 



Sf- 
T* = 
B = 



2 ms 

= i/r=io/zs 

10 nT 



n 0 = 10 piT Hz" 1/2 (43 nT Hz 
isotopically pure crystals 40 ) 

rj = 

530 fps 



1/2 for 



The total acquisition/integration time for each data point is the product of the binning 
time, At = 24.2 fis, and the number of cycles, N p . Hence we may plot the minimum 
detectable field, dB, against the integration time, T = N p At, as shown in Fig. 4(c). The 
sensitivity is defined by rj = 5B\ff and is plotted in Fig. 4(d), which is effectively 
constant. 
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